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In this talk I describe work, mostly done in collaboration with W. Simon some 20 
years ago [1,2], on multipole moments of static spacetimes. My purpose is to make this 
work, which lies at the interface of classical potential theory, conformal geometry and 
general relativity, known to mathematicians and to perhaps motivate them to have a 
look at the open problems which still remain. 

Our basic object will^be an asymptotically flat, smooth or analytic Riemannian 3- 
manifold (M,^) where M is diffeomorphic to R 3 \ B Ro with B Ro the Euclidean ball of 
radius R . In the chart given by this diffeomorphism the metric is required to satisfy 
the fall-off conditions 

9ij - Sij = 0(l/r), r 2 = 5 ij x t x J 

d~ 9ij = 0(l/r 2 ) (1) 
d% 3 = 0(l/r 3 ). 

It is often convenient and sometimes vital to have a stronger notion of asymptotic flatness 
which is defined via "conformal compactification" (CC) and goes as follows: There should 
exist a positive function Q on M, such that the metric = Q 2 gij extends to a smooth, 
or perhaps analytic metric on M = M U {r = oo} = MU {A}, where 

Q\a = 0, fti\ A = 0, (n ij -2g ij )\ A = 0, (2) 

with denoting covariant derivatives. There is a gauge freedom associated with this 
definition. Namely, if Q gives a CC, then Q = ufl, with cu\\ — 1, gives another one. Note 
that, with the above notion of asymptotic flatness, there exists a CC with a rescaled 
( "unphysical" ) metric being C° but not C 1 in general. Namely, pick VL = l/r 2 and 

x'^x'/r 2 ( "Kelvin inversion" ) (3) 

as coordinates on M. Then 

(4) 

with n H = x H /r' 2 , v! i = S^n'^ . Hence 

9'ij = + ($ k i - n' k n'J(6 t j - v! l v! 3 ){g ki - 5 M ) 

= 6 ij + 0(r'). (5) 

Clearly, if g^ is flat, gij is also flat, so there is a smooth, even analytic, CC in this case. In 
some applications where a smooth CC exists, the conformal factor Q can not be chosen 
to be smooth. As an important example take the spatial Schwarzschild metric. This is 
given by 

gijdx'dx 3 = ( 1 + — J dijdx'dx 3 (m > 0). (6) 

Then, taking x' 1 = x l /r 2 and Q = r /2 (l + mr'/2) 2 , we have that g^ = 5^-, but Q is only 
C 2 . Clearly, if we had started from 5^ rather than g^ as the basic metric this problem 
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would not have arisen. Fortunately, if is the spatial metric of an asymptotically flat 
static spacetime satisfying the Einstein vacuum equations, as is the Schwarzschild metric, 
there is a general way of writing Q = f'Q' so that Q' is smooth. An asymptotically flat 
vacuum spacetime, for the present purposes, is the manifold (R x M, ds 2 ) with ds 2 given 
by 

ds 2 = -e 2U dt 2 + g^dx?, (7) 
with U a function on M satisfying 

U = 0(l/r), dU = 0(l/r 2 ), d 2 U = 0{l/r 3 ). (8) 

We now rewrite as 

9H = e~ 2 %- (9) 
In the case of the Schwarzschild metric we have that 



e 2U = 



/o \ 4 / 2 \ 2 

m/2r\ „ / m 



l + m/2rj ' , J \ 4r 2 



7« = 1 - 7-2) 5 ^ r > m / 2 - ( 10 ) 



(Thus Rq has to be taken > m/2.) It easily follows that 7^ has a smooth CC with 
smooth conformal factor. The Schwarzschild solution is a solution of the Einstein vacuum 
equations on R x M. These equations are equivalent to 

A~U = 0, (11) 

H ij [T\=2U i U j , (12) 

where A~ is the Laplace Beltrami operator on (M, 7^) and ^ [7] the Ricci tensor of 
7ij. The equations (11,12) should be viewed as the "source-free static Einstein equations 
on an exterior domain". In particular, (U, 7^) could have smooth extensions to R 3 so 
that ds 2 satisfies the static Einstein equations on R x R 3 with an energy momentum 
tensor corresponding to some "reasonable" matter model. If, for example, this matter 
model corresponds to a perfect fluid with density p, pressure p and an equation of state 
of the form p = p(p) satisfying a certain differential inequality, it is known [3,4] that the 
spacetime is spherically symmetric (whence Schwarzschild in R x M), and thus there is 
just a one-parameter family of solutions, parametrized, say, by the central density. In the 
trivial case that the spacetime is vacuum throughout R x R 3 it is known [5] that ds 2 is 
the Minkowski metric. (If the matter model describes an elastic solid such as the earth, 
the solutions will in general be neither Schwarzschild nor flat space, of course.) The point 
of restricting ourselves to an exterior domain is that we do not care about what sort of 
matter may be "inside". Then there is a larger class of solutions. How large, in fact, this 
class is and how it can be conveniently characterized, is precisely the question we are 
asking. On physical grounds one expects that there should be the same freedom as in the 
Newtonian theory, to which we now turn. The Equ.'s (11,12) are now replaced by 

A~U = 0, U = 0(l/r) (13) 
72y[7] = 0. (14) 
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Thus, since dim M — 3, the Riemann tensor of 7^ has to vanish. By the simple connect- 
edness of M it follows that 7^ = 5ij in suitable coordinates, and thus there remains 

A s U = 0, (15) 

that is to say the Laplace equation with respect to the Euclidean metric. It is convenient 
to replace (15) by 

AsU = -4np, (16) 
where p G C£°(R 3 ), whence U has to be of the form 

U(x) = I -^-dx'. (17) 

J-R? \X - X'\ 

It follows [6] that U has an expansion of the form 

~ _ M Mix' 1 M ij x i x j M h ik x h . . . ./•" 

U ~ — + ~^3~ + 2r 5 + • • • + k\r 2k+1 + "' ( ' 

where the constants M il ^ ik are symmetric and tracefree, which converges uniformly in 
M (see [6]). As a warm-up for general relativity we outline an independent argument for 
this classical result, as follows: One first shows that 

~ k M- t* 1 ■r ii 

"-;1 m +od/o (w) 

for some sufficiently large k, say k = 2, and that this relation can be differentiated at 
least twice (i.e. any partial derivative of the left-hand side is 0(l/r k+2 ), a.s.o.). Then 
pick the standard CC. Now define an unphysical potential U by U — VL~ l l 2 U . This has 
the form 

U = M + + ^Mijx'W + 0(r' 3 ). (20) 

Thus U is C 2 ' a (0 < a < 1) in B^m — M. Next observe that U again satisfies the 
Laplace equation. This is seen by first recalling the conformal behaviour of the conformal 
Laplacian, i.e. 

a 7 - 1) u = ir 5 / 2 (a~ - ^ju, u = n- 1 ' 2 u (21) 

and specializing to the case where both 7^ and 7^ are flat. We now have that U is a 
C 2,a -solution of 

A S U = in M. (22) 

But it is well known that solutions to the Laplace equation have to be analytic [7] . Thus 
U has a convergent Taylor expansion at the origin, which ends the proof. 

Since analyticity, by [8], is a general property of solutions to elliptic systems with 
analytic coefficients, one can try to use a similar line of thought for the relativistic case. 
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The main obstacle is that the equations (11,12), far from being invariant under the 
conformal rescaling 

u = n-^u, ll3 = n%„ (23) 

become formally singular at the point A where Q vanishes. Before showing how this ob- 
stacle can be overcome we have to perform the analogue of the first step in the Newtonian 
case, namely the analysis of the first few terms of the expansion of the quantities (U, 7^) 
in powers of 1/r. This has been done, up to arbitrary orders in 1/r, in the sequence of 
works [9,10,1]. One could obtain sufficiently detailed information about the general term 
in order to show: For every k there exist coordinates x % so that, after the standard CC, 
7ij is smooth. For example, for k — 3, one finds 

7, = % - V - - ^P 1 + ^^(-V + + 0(1//), (24) 

& = ^ + ^ + i^ + ^ + (l/^) (25) 

for some constants M, Mj, Mjj (Mjj symmetric and tracefree) and where Xi = SijX J ' . 
These equations show in particular that, at least at this order, the relativistic solution 
has no more free parameters than the Newtonian one, and this remains true for arbitrary 
high orders. 

Very loosely speaking, whereas there are harmonic functions solving Equ. (11) with 
the appropriate boundary conditions which are parametrized by multipole moments, the 
metric 7 is essentially determined by the r.h. side of Equ. (12). This, in turn, is related to 
the fact that IZij determines the curvature tensor in 3 dimensions and that zero curvature 
implies a flat metric. Within the iterative scheme by which the system (11,12) is solved 
this means that, for each order in 1/r, "harmonic" contributions to 7^ have to be pure 
gauge. 

Armed with this information one can now study a suitable "unphysical" version of 
Equ.'s (11,12). It is not clear how to use the standard CC for this purpose. The basic 
observation is that the expansion of U, as in (25), also says that U 2 is arbitrarily smooth 
in Kelvin inverted coordinates and that, when M 7^ 0, (fi, 7^) satisfies conditions (2) 
with = U 2 /M 2 . It is instructive to first see the effect of this CC in the Newtonian case. 
In that case we have that 1Z = so that AU = and U = M. Now Equ. (21) implies 

TZ = 0. (26) 

Conversely we find from (A — TZ/8) - const = that, outside A, A(f2 _1 ) = 0. Equivalently 

QAtt = (27) 

Of course, the unphysical metric 7^ is now not flat any longer in general. It satisfies 

- nilij = - ^gijAQj . (28) 
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The relation (27), as it stands, does not give rise to a regular elliptic equation for Q since 

° - 2— < 29 > 

is formally singular at A. Using (2.8) we find that 

<t j = -SKijW (30) 

which further implies 

Aa = milijTCK (31) 

Now (31), together with Af2 = a, does form a regular elliptic system. Using the theorem 
of Morrey [8], it would now follow that (Q, cr) are analytic provided 7^ is. This, in turn, 
follows from the following differential geometric 

Lemma: A conformally flat metric with zero scalar curvature is analytic. 

Proof: By 1Z = and conformal flatness there holds 

DkR-ij = DilZkj- (32) 
Taking D k of Equ. (32), commuting derivatives and using 

T^kju = ^i[k9j]e — ^e[k9j]i (33) 

and the Bianchi identity we obtain 

AKij = mfKjt + gijlZull™. (34) 

Now recalling that the Ricci tensor of a metric gij, when written in harmonic coordinates, 
gives an elliptic operator for 7^, we see that 

together with (34) yields an elliptic system for (7^,0"^) with analytic coefficients. Thus, 
by [8], the Lemma is proved. 

In the case at hand we have that 1Z vanishes and, furthermore, that 7^ is conformally 
flat. In fact, from (29) we can directly obtain 

TZ im = 0. (36) 

Combining the analyticity of 7^ with the system (30,31) we obtain analyticity of the pair 

The relativistic situation, given by Equ.'s (13,14), is similar, but more complicated. 
Instead of fi = U 2 /M 2 it is more convenient to take 

4 U 

Q = — sinh 2 - (37) 
M 2 2 v ; 
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as conformal factor. As a field variable it is useful to take, instead of f2, the quantity p, 
defined by 

p = tanh 2 — . (37') 
Then, again, it turns out that 1Z — and 



a 3 PiP 

Ap = a:=- ' ( 38 ) 

2 p 



whereas (28) gets replaced by 



-p(l - p)Hij =Pij - ^YijAp. (39) 

Clearly the metric 'jij can not any longer be conformally flat. In fact, taking the "curl" 
of (39), Equ. (36) of the Newtonian theory gets replaced by 

(1 - p) Ui[j; k ] = 2TZi\jp k] + ^i\jTZ k]e p e . (40) 

As a side remark we add, that, from the above equation it is not hard to show that is 
conformally flat iff 7^) comes from the Schwarzschild solution. 
For a we obtain, instead of (31), the relation 

Aa = 3p (1 - pf + ZKijpy. (41) 

Apropos Equ. (41) we wish to point out that, with the substitution 

and using (40), there results the identity 

AK = ^K 5 + ^K- 7 n im K^ k \ (43) 

which plays an important role in black hole uniqueness theory (see e.g. [11]). 

In any case, we can now consider the system consisting of (35), (38), (41) and the 
equation of the form 

Aa y = . . . (44) 

obtained by multiplying (40) by (1 — p) _1 (which is non-zero near A) and taking D k . 
This furnishes an elliptic system with analytic coefficients for the variables (p, 7^, a, a^) 
(that-is-to-say if harmonic coordinates are used in Equ. (35)). We conclude that (^,7^) 
is analytic. 

In the Newtonian case, the multipole moments of the solution were nothing but the 
Taylor coefficients at A of the unphysical potential U in the standard CC. In particular 
this implies that these moments determine the physical solution uniquely (which however 
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already follows from the expansion (18)). In the relativistic case, and using the CC above 
based on Q = U 2 /M 2 , it is not immediately clear whether the moments determine the 
solution. In fact, the expansions in (24,25) use a specific coordinate condition. In the 
unphysical picture one would like a "gauge invariant" definition of multipole moments. 
Luckily, Geroch in [12] came up with such a definition. It goes as follows: Consider the 
following recursively defined set of tensor fields built from (U, Sl,~fij) 



p 


= u 

— TT. 






P 






is + l 


= C 


s(2s-l) 
D P 1Z P 

-^ts+l 1 11— Is 2 /V lll2- 1 13---1 


3 + 1 



(45) 



where C denotes the operation of taking the tracefree symmetric part. Now define 



^ii...i s ■ Pii...i s Ia ' 



(46) 



Clearly, in the Newtonian case and using the standard CC for which IZij = 0, this 
definition coincides with the standard one. As for the gauge dependence, when 



n = un, u = uj- 1/2 u, 

with u\ A — 1 it turns out [13] that 



(47) 



Mi 




(2s - 1) . . . (2r + l)(-iy- s M tl ... ia b tr+ i ...b is , (48) 



where bi := \uJi . Note that only first derivatives of uj at A enter the transformation law 
(48). Furthermore the dependence of M^...^ on bi is exactly the same as the change of the 
Newtonian moments under a translation of the physical Euclidean space by the vector 
bi. We can now come back to the conformal gauge with Q = U 2 /M 2 . Here it follows 
from (45) that U = M so that the moments are equivalent to M, together with C of the 
derivatives of IZ^ at A. But it is not difficult to infer from the above equations that these 
data determine (Q, 7^) uniquely whence (U, 7^) is also determined. We remark, finally, 
that the M^. ^-terms in the physical expansion, which we have written out for k = 2 in 
(24,25), coincide with the Geroch moments from above. 

There are two problems which the above analysis leaves open. The first one is to give 
an analyticity proof in the case where M is zero. Although the "physical" expansion of 
(U, 7ij) seems to be completely insensitive to whether M is nonzero or not, the unphysical 
situation, because of the choice of conformal factor makes vital use of the nonvanishing 
of M. The second open issue is that of convergence. Namely, given a set of moments 
Mil. ..i 3 > s — 0,1, . . .: what is the condition on the behaviour for large s so that a static 
vacuum solution having these as multipole moments exists? The Newtonian situation is 
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clear: the moments have to be such that the Taylor series having them as coefficients 
converges. In the relativistic case, the Taylor coefficients at A of the quantities (fi, 7^), 
say in Riemannian normal coordinates, depend on these moments in a nonlinear fashion. 
In particular it is not even obvious whether a solution exists for which only finitely many 
moments are nonzero. 
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